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■ Abstract 



The standard generators of tridiagonal algebras, recently introduced by Terwilliger, are shown to generate 
a new (in)finite family of mutually commuting operators which extends the Dolan-Grady construction. The 



> 

involution property relies on the tridiagonal algebraic structure associated with a deformation parameter q. Rep- 
resentations are shown to be generated from a class of quadratic algebras, namely the reflection equations. The 
spectral problem is briefly discussed. Finally, related massive quantum integrable models are shown to be super- 
integrable. 
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1 Introduction 

O ' 



In the Liouville- Arnold sense, a classical Hamiltonian system with N degrees of freedom is called completely 
integrable if there exists a set of N functionally independent globally defined integrals of motion in involution, 
i.e. mutually commuting regarding to the Poisson bracket. It is said to be superintegrable if it admits more than 
N integrals of motion. By analogy, a quantum system with N degrees of freedom is completely integrable if there 
exists a set of A^ — 1 operators, say TL n , n = 1, N — 1, together with the Hamiltonian TL which are mutually 
£h commuting. If there exists p— additional operators X m , m = 1, ...,p where 1 < p < N such that 



[H n ,n m }=0 and [H,I m }=0 (1) 



for all n, m, this quantum system is said to be superintegrable. Although not necessary for superintegrability, the 
additional operators T m can be in involution too i.e. [I n ,I m ] = 0. 

The purpose of this paper is to introduce a new (in)finite set of mutually commuting quantities, whose 
involution property will be indentified with the defining relations for the tridiagonal algebras recently introduced 
and studied by Terwilliger £Q (see also We will also show that the underlying structure is closely related 

with U q i/2 (sfa), and coincides with the Dolan-Grady one 0] for 5 = 1. As a consequence, we will show that related 
quantum integrable models (open spin chains, sine-Gordon field theory,...) enjoy superintegrability. 

This paper is organized as follows. In Section 2, after recalling some definitions we propose a new realization of 
tridiagonal pairs in terms of 1/^1/2(5/2) generators. In particular, the defining relations of the tridiagonal algebra 
are shown to be invariant under the coproduct homomorphism of U q i/2{sl2)- Based on the tridiagonal algebraic 
structure, we propose a finite set of conserved quantities T n in involution. In Section 3, an alternative construction 
of the conserved quantities I n is described. Their generating function is expressed in terms of general solutions 
of a class of quadratic algebras, namely the reflection equations. In Section 4, , we argue that the corresponding 
spectral problem is related with a system of N partial q— difference equations. For iV = 1, the eigenvalue is 
expressed in terms of solutions of Bethe equations and enjoys a remarkable symmetry property. Examples of 
massive quantum integrable models with (in)finite degrees of freedom are considered in Section 5, where the 
quantities I n are explicitly identified. Concluding remarks follow in the last Section. 

1 baseilha@phys.univ-tours.fr 
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2 Tridiagonal algebraic structure and involution 



Tridiagonal algebras have been introduced and studied in [UEID]) where they first appeared in the context of P— 
and Q— polynomial association schemes. Wc will here consider a slightly simpler form of the tridiagonal algebraic 
structure, which nevertheless possesses all interesting features for further analysis. For more informations about 
these algebras, we report the reader to £Q where precise definitions can be found. Let us consider the tridiagonal 
(associative) algebra T with unity generated by two operators (called standard generators) acting on a vector 
space V, say A : V — ► V and A* : V — + V, subjects to the tridiagonal relations 

[A, A 2 A* + A* A 2 — (q + q~ 1 )AA*A - pA*] = 0, 

[A*,A* 2 A + AA* 2 - {q + q- 1 )A*AA* - pA] = (2) 

where q is a deformation parameter (assumed to be not a root of unity) and p is a fixed scalar. According to 
[PP, Theorem 3.10] A, A* is a tridiagonal (TD) pair as defined in Definition 1.1], which complete classification 
is still an open problem. Among the known examples of TD pairs, one finds a subset such that A, A* have 
eigenspaces of dimension one. These are called Leonard pairs, classified in 0. In particular, they satisfy (for 
details, see |B]) the so-called Askey- Wilson (AW) relations first introduced by Zhedanov 0. Other examples of 
TD pairs can be found in 0|S]: for p = in which case reduce to q— Serre relations or for q = 1, p = 16 which 
leads to Dolan-Grady relations @]. 

Here, our aim is to construct a new family of TD pairs in terms of £7^1/2 (s^) satisfying J5J. To do that we 
introduce five generators Q±, Q ± and H subjects to the algebraic relations 

q- 1/2 Q±Q± - q 1/2 Q±Q± = o , 

_ _ „±2H _ i 



=F * ^ T ^± q l/2_ q -l/2 



q eH Q±=q ±e Q±q eH , q eH Q± = q ±e Q±q eH (3) 



and the q— Serre relations 



Q 3 ± Q T - (1 + q + g- 1 )g 2 b g T Q± + (1 + q + q- 1 )Q ± Q T Q 2 ± - Q T Q± = , 
Q 3 ± Q T - (1 + q + q- 1 )Q 2 ± Q T Q ± + (1 + q + q- l )Q ± Q T Q 2 ± - Q T Q± = . 

Using the standard definitions it is an exercise to show that Q±, Q± and H actually generate the quantum 
affine Kac-Moody algebra U q \/i{sl2). Based on the results of JHIi we propose 2 the following TD pair: 



together with the scalar 



A = -Q+ + Q_ and A* = g_ + — Q, , (4) 

co c 

(g 1/2 + T 1/2 ) 2 

p = . (5) 

co 

It is straightforward to check using the defining relations above 10 that Q satisfy the tridiagonal relations . 
Furthermore, we have found that the slightly more general TD pair A = A + e + q H , A = A* +e-q~ H also satisfies 
(J2J for any values of the parameters e± . This is not surprising as we will see later on. 

The Hopf algebraic structure of [/^^(s^) can now be used to construct a whole family of TD pairs. Indeed, 
the coproduct A : U q i/i(sl2) — > E/ ? i/2 (s^) x U q i/2(sl?) associated with Q is such that 

A(Q±) = Q±®I + q ±H ®Q±, 
A(Q ± ) = Q ± ®I + q^ H ®Q ± , 

Afef) = q H ®q H . (6) 

2 This realization is different from the one proposed by Terwilliger and Ito in |S| in which case p = 0, i.e. reduce to q— Serre 
relations. 
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More generally, one defines the N— coproduct : 17 1/2(3/3) — > U q i/2{sh) ® • • • <8 U q i/2{sl2) as 

AW = (id x • • • x idx A) o A^- 1 ) (7) 

for TV > 3 with A< 2 > = A, A* 1 ) = id. The opposite iV-coproduct A'W is similarly defined with A' = a o A 
where the permutation map cr(a; ®y) = y®x for all aj, y 6 C/ ? i/2 (sZa) is used. Due to the homomorphism property 
of the coproduct (or its opposite) A(xy) — A(x)A(y) the algebraic relations © are invariant under the action of 
A (or A'). Consequently, the family of TD pairs (acting on tensor product of irreducible representations) defined 
as 

AW=AW(A) and A * W = A^(A*) for N > 1 (8) 

with @ satisfies the tridiagonal relations 3 (J2J with JSJ. It is interesting to notice that the TD pair A, A can 
also be generalized as in ©. The corresponding tridiagonal algebra TW is a (left) coideal subalgebra of 
U q \/z{sl2). Indeed, one has A(a^) € £^1/2^2) <8> for all a £ {J, A, A }. As a result, this class of coideal 
subalgebras closes under the tridiagonal algebraic relations 

An interesting problem is whether the tridiagonal algebraic relations J5J might provide the condition of exis- 
tence of a finite set of mutually commuting quantities. Having in mind the construction of Dolan and Grady g] 
associated with (j2J for q = 1, whose generators can be expressed in terms of the loop algebra of s/2 |12l I13| . this 
problem seems rather natural and solvable. Inspired by @], JHI we propose the following finite set of quantities 
in involution, denoted T^n-i witn n=l,...,N, generated by the TD pair A W , A* {N) ® with (gj: 

xW = kA w + «*a* w , 

Zf) = K ([[ A W,A* m ] 9 ,AW], + pA ,(Ar) )+ fi *([[A* W ,A W ],,A ,W ],+/'A (w) ), 

4 W) = 9 9 (9) 

where the g— commutator [X, Y] q = q x l 2 XY — q~ x l 2 YX has been introduced and k, n* are arbitrary parameters. 
As shown in the first part of this Section, the TD pair 

A (jv) A *W 

satisfies the tridiagonal algebraic relations © . 
Then, it is an exercise to check, for instance, that [x[ N ,T^ N ] = for any N. Note that the special case n = n* 
was considered in [TJ]|. For general values of N, it is not difficult by analogy with 0] to guess the explicit form 
of the N independent quantities for n = 3,...,N (higher order polynomials in the TD pair) which are mutually 
commuting. However, due the q— commutator the proof of mutual commutativity of the corresponding quantities 
becomes quickly rather complicated. Consequently, in next Section we propose an alternative method to derive 
these quantities. 

To see the truncation of the hierarchy at the order N let us focus on the smallest values of N. For N = 1, using 
the evaluation homomorphism n v : U q i/2(sl2) — ► 1/^1/2(3/2) (i> 6 C is sometimes called the spectral parameter) 

tt v [Q±\ = vS±q ±S3/2 and tt v [Q ± ] = v^S^^ 2 7r„ [H] = s 3 , (10) 

where {S±, S3} are the fundamental generators of the quantum algebra U q i/2 (sh) satisfying the relations [S3, S±] = 
±S± and [S+, S-] = (<f 3 - q- S3 )/(q^ 2 - q~ 1/2 ), one finds that A (1) , A* (1) satisfy the AW relations 4 

A ,(1) >+aWa' (1)J -(? + 9 - 1 )A ,{1, aWa ,(1) = P A«+ W A* (1) (11) 

with (JSJ and 

UJ = -{v 2 q- 1 ' 2 + v- 2 q l ' 2 )w^/c . (12) 

3 Note that an other family of TD pairs can similarly be obtained using instead the opposite coproduct A'. 
4 Here we use the notations of |§] for p = p* , 7 = 7* = 0, r? = rj* = 0. 
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Here, w^' = {qi +1 / 2 + q J x / 2 ) denotes the eigenvalue of the Casimir operator of [^1/2(3/2) in the spin— j 
representation. In other words, the TD pair A*- 1 -*, A*'" 1 '' is a Leonard pair (see [5] for details). Consequently, one 
can not construct an independent quantity trilincar in A (1) , A* (1) in the spirit of 0], due to Ijlll) . Appart from the 
Casimir operator of the AW algebra y\ (quartic in the Leonard pair), there are no other independent quantities 
commuting with l[ directly built from the TD pair. 

For N = 2, we have checked explicitly that the TD pair A*- 2 -*, A**- 2 "* does not satisfy the AW relations As 
before, we expect that an independent quantity 1^ of fifth order in f^ 2 \ A**" 2 "* does not exist. As will be shown 
in the next section, this fact relies on the existence of algebraic relations of fifth order in the TD pair. More 
generally, there exists a set of algebraic relations generalizing the AW ones responsible of the truncation of this 
integrable hierarchy. In the next Section we will propose a general procedure to derive all l^n—i f° r 77 = 1' •••> 
all algebraic relations generalizing I jlljl and the involution property of the charges I^n-i- ^ * s based on the 
generalized quantum inverse scattering approach. 



3 Generating function 

In the spirit of the quantum inverse scattering method, we would like to construct an object which would provide 
the generating function of the above mentioned mutually commuting quantities ^n-i! n = 1, N. Following the 
results of let us consider the following quadratic algebra (called reflection equation) which was introduced 
by Cherednik [T4"] : 

R(u/v) [K{u) <g> 1) R{uv) (J (8> K{v)) = (M <g> K(v)) R{uv) (K(u) ® I) R(u/v) . (13) 

This equation arises, for instance, in the context of quantum integrable systems with boundaries |15j . We report 
the reader to the literature on the subject for more details. For our purpose, we restrict our attention to the 
trigonometric R— matrix R(u) which solves the Yang-Baxter equation. In the spin— | representation of U q i/2(sl2), 
it reads 

R(u) = ^2 Vij(u) o-i® <jj , (14) 
ije{o,3,±} 

where 

oj 00 (u) = 1(^/2 + i)( u -g-Va u -i )) W33 (u) = i(« 1 /2_i)( u + g -Va u -i )) 

= w-+(«) - <z 1/2 - <T 1/2 

and tjj are Pauli matrices, o~± — (a~i ± ia^jl. As shown in |15| . the solution of is not unique: a family of 
solutions of the reflection equation can be obtained from an "initial" solution, say K^°'(u), using the so-called 
"dressing" procedure. Indeed, given the Lax operator L(u) which satisfies the Yang-Baxter algebra 

R(u/v)(L(u) <g> L(v)) = (L(v) ® L{u))R{u/v) (15) 

then one obtains a family of "dressed" reflection matrix for any parameter k G C 

KW(u) = %(ufc) • • • L^u^K^i^L^uk" 1 ) ■ ■ ■ L^uk- 1 ) (16) 

acting on the quantum space = ® j = l which also solves Following the analysis of jTH], let us then 

introduce the elementary solution of the "dual" reflection equation 5 given by 

/ uq l / 2 K + u- l q- l / 2 K* \ , 17 , 



5 This equation is obtained from ll.il by changing u — ► u 1 , v — * v 1 and K(u) in its transpose. 
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where k* = k 1 for any k£C. Together with JTBJl, it constitutes a basic element for the construction of mutually 
commuting quantities. Indeed, for any values of the spectral parameters u,v 

[t w (u),t^ N \v)] = for t^ N \u)=tr {K + (u)K^(u)} (18) 

where tro denotes the trace over the two-dimensional auxiliary space. We refer the reader to ^1] f° r details. We 
are going to show that, given N and the choice K^(u) = (a + /c + cr_)/(g 1 / 2 — q^ 1 ^ 2 ), the set of N mutually 
commuting quantities ^n— 1> n = 1, N described in previous Section are generated from . In the following, 
we introduce a second spectral parameter v — kq 1 ^ 4 and we take the fundamental solution of the Yang-Baxter 
algebra (L— operator) in terms of U q i/?{sl2) generators: 



L ( U )-[ q \ q ^- q -^)S + uqiq- s ^ 2 -u- l q-^q s ^ 2 ) ' 



ss/2 _ u -l q -\ q -s 3 /2 ((? l/2 _ q-l/2) S _ 



Let us assume 



KW(u)= cr j ®nf ) (u) (20) 

je{o,3,±} 

where Qj(u) G Tun[u\ fiS N \ A*^- 1 ) are Laurent polynomials of degree — 2N < d < 2N in the spectral parameter 

u, the TD pair being subject to the tridiagonal relations J2J- For N = 1, using the evaluation 

homomorphism the result for K^ 1 ' (u) coincides with the - called "dynamical" - solution of the reflection 
equation proposed in ^H] (see also JHj)- Explicitly, in terms of the Leonard pair A^, A*^ it is given by l(20|l 
with 

nW^+n^w = V/2 A W- u -i r V2 A *(i) t 
^\u)-n^(u) = V/ 2 a* (1) - h-V 1/2 a (1) , 

= ql/ 2 u 2 + q-V 2 u- 2 , [A* (1) ,A ( \ , UJ 



c (q 1/2 -q- 1/2 ) q 1 ' 2 + q- 1 ' 2 (q-q- 1 )' 
(1) _ g l/2 u 2 +g -l/2 u -2 Co[A W,A*(\ cpc 

W - ( 9 l/2- 9 -l/2) + gl/2+ g -l/2 + ( (Z - 9 -l)- 

^From this expression, we deduce immediately t^ N=1 ^ — {qu 2 — q~ l u~ 2 )T[ with 

= K A (1) +k*A* {1) , (22) 

where the evaluation homomorphism 1)10(1 has been used. In particular, I± is linear in the generators A^, A* , 
as anticipated in the previous Section. There are no higher quantities X^n-i for n > 2 generated in the expansion, 
a phenomena which relies on the existence of the AW relations (|llfl which follow from the reflection equation |10j . 

For N — 2, the quantities ,2^ can also be exhibited in the expansion of K^ 2 \u). To show that, we first 
use eqs. J2J, © and (JBJ) to derive 



[A ( 2)) [A^A^J^ = [AW, [AW.A'W] J rl ® J + «*® [A«, [A« A*«] J ? _ 

(1 ~ g ' 1) Q + + (l-g)Q UAWA*« 
co 

ii^Q + + (l- g - 1 )Q )®A*WAW 

c / 



CO 



qQ+Q^_ - q- 1 Q^Q+)q- H ® A* (1) 



" + Q+Q-)g g + (1 g 1} (Q-Q+ + Q-Q+)g g ) ® A (1) , (23) 



Co c 
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and similarly for [A* (2) , [A* (2) , A (2) ] ] _ lf exchanging q -> q' 1 , A (1) A* (1) , Q+ <-> Q + and Q_ <-> Q_ in the 



above expression. Using the evaluation homomorphism l|l(J|) . we replace the AW relations IjllH in the two first 
terms of <|23|) . A straightforward calculation finally shows that the diagonal entries of K^ 2 \u) can be written in 
terms of <|23[) . We obtain for instance 



n { 2 \u) + n^\u) = (qu 3 + uq l ' 2 f{v) + u-^A^ - (q^u- 3 + u^q- 1 ' 2 f{v) + u)A* {2) 

/ ^1/2 _i_ „-l/2\2 

— 1/2 C ° ^ (-[A^A^A^U-* + ^ c } A- 



(2) 



where 



2(z,v^ + t ,-y/^ w o) 

= g V2 + g -i/2 ■ ( 24 ) 

The expression for fi^^u) — D,^\u) is derived similarly, and corresponds to the substitution A^ 2 -* «-> A in the 
above equation. The mutually commuting quantities can now be read off from the coefficients of the expansion 
in the spectral parameter u of 11811 . Taking the trace over the auxiliary space of (1201) for N — 2 together with the 

(2) 

expression for S2q , we find 

t^(u) = (g 3 /V - q- 3 ^ + Ctf,u-*;v))7?> + C ° ( g T/ 2 2 ; 9 g -!/ 2 ' 2) 4 2) , (25) 

where l[ 2 \l^ coincide exactly with @ for N = 2 and the Laurent polynomial C{u, v ) directly follows from 
H24|) . The property l|18|) implies that the two quantities l[ 2 \x^ are mutually commuting. This provides an 
alternative derivation to the one proposed in previous Section which is solely based on the fact that A*- 2 -*, A**- 2 " 1 
satisfy the tridiagonal algebraic relations |J2}. 

This technical procedure can be clearly applied for general values of N. In this case, it is important to notice 
that the leading terms in the asymptotic u — > oo yields to 

U 2-2N q (l-N)/2 K (N) (u)l J * ^ C ° (9l/2 T 1/2) ^ (M) 



Similarly, the asymptotic u — ► — oo gives an analogous expression exchanging A^ «-> A . We immediately 
recognize the TD pair defined in (JSJ) which obeys J2Jl as shown in previous Section. Then, the mutually commuting 
quantities Q are obtained from the expansion of t^ N ' (u) in the spectral parameter u which can be formally written 

N 

^)( M ) = ^CW 1 ( U 2 "; M - 2 ",.., W 2 , U - 2 ;«)IW 1 , (27) 

71=1 

where C^n-i are Laurent polynomials in u, w. According to 1)180. they are in involution i.e. [^2n-i'-^2m-i] = ^ 
for any n, m > 1. Such property is actually encoded in the tridiagonal relations J2J. It should be stressed that 
the above approach does not depend on the choice of a representation for (|10|1 . Consequently, finite, infinite 
dimensional or cyclic representations of TD pairs and related quantities like T^n—i can ^ e easily obtained. 



4 Related spectral problem 

In the context of quantum integrable systems (see next Section), an important problem is the diagonalization of 
lfTT|) which leads to study the spectral problem associated with © that we shall briefly discuss now. We write it 
as 

4V-i*l N \zu---,ZN)=AW 1 *i N Xz 1 ,...,z N ) , for n=l,...,N, (28) 
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where (zi, zn) denote the eigenfunctions. For q not a root of unity and general values N, we expect 
that this problem can be "algebraized" 6 i.e. part of the spectrum corresponds to multivariable polynomial 
eigenfunctions. Here we do not attempt to study in details the corresponding system of N partial (/—difference 
equations, which goes beyond the scope of this paper. Instead, let us focus on the case N = 1 which can be solved 
easily following (171 IT] (see also JOI)- Expressed in the weight basis, the generators of J7 g i/2(s^) leave invariant 
the linear space (of 2j + 1 dimension) of polynomials F(z) of degree 2j. The lowest/highest weights are identified 
with Fq — 1 and i*^ = z 2j , respectively. One has 

S + F(z) = ~ (gl/2 _* g -i/2) (l- j F{<f' 2 z) - cfFiq-^z)) , 

= z{ql/2 [ q _ 1/2) (F( q V*z)-F{q-V*z)) . 

In this basis of single- variable polynomials, the spectral problem for © with N = 1 reduces to a single second- 
order q— difference equation of the form 

a (1 \z;v)^(qz) + d {1) (z;v)^(q- 1 z)-v^(z;v)^(z) = A^ } *(z) . (29) 

Indeed, we should keep in mind that - at this special value of N - all the higher charges can be written in terms 
of the first one due to the AW relations Ijlljl. In (|29|l . the coefficients are given by 

1,1 CQKV~ 1 q~^ 2 z~ 1 — nvq~ 3 ^ 2 z 



oW(z;«) 
d {1) {z;v 



K*v~ 1 q 3j / 2 z - c K*vq J / 2 z~ 

cote 1 / 2 - <r 1/2 ) 



m . . (c kv x q ] / 2 - c K*vq^ 2 )z 1 + (k*v x q J / 2 — KVq^ 2 ) 
V {Z;V) = C^-q-^ 2 ) 



(30) 



and 



M 



*{z) =Y[(z-z m ), (31) 



m— 1 



where z m denote the roots of the polynomial. Dividing 1(29(1 by ijSH , the l.h.s. of 1)29(1 gives a meromorphic 
function of z and the r.h.s. becomes a constant. Singularities of the l.h.s must be cancelled. They are located at 
z = 0, oo and z = z m . For z = 0, 1)29(1 vanishes identically. For z = oo and general values of /c, k* on gets the 
constraint M = 2j. For z = z m , one finds the following system of Bethe equations: 



M=2j 



p\= n q —^^ *» i =i>~>v pa) 

to which corresponds exactly 2j + 1 polynomial eigenfunctions. Comparing the constant terms of both sides in 
<(29|) one obtains 

Af 3 = - V j/2+1/2 + wu^/a-i/a) ^ Zm / Co . (33) 

m— 1 

Let us also mention an other way to solve the spectral problem ((28(1 for TV = 1, which exhibits an interesting 
link with Askey- Wilson q— orthogonal polynomials. These Laurent polynomials are symmetric in the variable y 
and defined by 

p-M-Ar^^-j^i.,,) (34) 



e We use the definition of |17l . 
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where 4$3 denotes the basic q— hypergeometric function, n is an integer and a, 6, c, d are arbitrary parameters. In 
particular, the zeros of the AW polynomials are determined by the system of Bethe equations [*T7] 

{Vk - a){Vk - b)(y k - c)(y k - d) = -pr - Vi){qykVi - 1) for fc = 1 n ^ 

[ay k - l)(by k - l)(cy k - l)(dy k - I) l= \fu k (Vk ~ m){VkVl ~ (?) ' 

As shown in (see also [J] for finite dimensional representations), the Leonard pair A (1) ,A* (1) can be written in 
this basis. Introducing the so-called Askey- Wilson second order q— difference operator 

D - i(y)(j - I) + av^Kr- 1 -/) + (! + abcdq- 1 )! (36) 

where r(y) — qy and 

p{ s_ (1 - ay)(l - by)(l - cy)(l - dy) 
^ (V> (l-y 2 )(l-qy 2 ) 

one can define the representation ir' such that 

n'[A^] ~(y + y-^Pniy) and tt'[A* (1) ] ~ DP n (y) . (37) 

Using the explicit form (|36|l in (|28(l . setting abed = q, cq = 1 and rescaling the parameters k, k* one obtains a 
second order q— difference equation of the form 

K*m?n(qy) + K^Ctf-'jPnfe -1 !/) - + CO/" 1 ) - ~<V + 2T 1 ) - 2))P n (y) = A^P^y) . (38) 

Similarly to the case discussed above, the eigenvalues Aj can be written in terms of the zeros yi^yf 1 of the AW 
polynomials determined by 

As both representations are available, let us consider the first (finite dimensional representation) 131|) . Then, 
it is worth important to notice that the Bethe equations l-Sl'l) as well as the eigenvalue A^ are invariant under 
the symmetry transformation 

q <-> q^ 1 , v <-> u _1 , k <-> k* . (39) 

Although we will not discuss this remarkable symmetry for general values of N, we expect it will be preserved for 
higher values of N too. 



5 Examples 

For several quantum completely integrable models, it is not difficult to construct non-local integrals of motion in 
terms of l^n-i- Below, we list some examples (see also [TU]1. 

• Integrable (XXZ) spin chain models: In the context of Sklyanin formalism, integrable boundary con- 
ditions for a system of finite size are encoded in the representations of the A"— matrix K^ N \u). From previous 
analysis, the problem is reduced to classify all possible (tensor product) representations of A^, A*^' associated 
with the tridiagonal algebra J5J). For q not a root of unity, the representations and (JSJ) are proposed. 

All known examples of integrable boundary conditions can be indeed seen as special cases of (jHJ. For instance, 
let us denote tt^ the spin— j irreps. of U q i/a(sl^) with 11U|) . For N = 1 and the trivial representation we can 

define 7r^=°)[A (1) ] = e_, 7r^=°) [A* (1) ] = e+. Then (id x 7T (i=0 ) ) [K W (u)] coincides with the non-diagonal solution 
of the reflection equation for the open XXZ spin chain with M sites and non-diagonal boundary conditions found 
in |18j . The method of deriving the "boundary quantum group" generators |19j proposed in |2()j was recently 
applied to this model [21] ■ It is then easy to see that the so-called "boundary quantum group" is actually the 
tridiagonal algebra (0), the boundary generators being nothing but a TD pair of the form A< M ), A W in the 
spin— g representation. According to Section 2, this model possesses M mutually conserved quantities given by 
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Q for k — k* = 1. Interestingly, the corresponding spectral problem gives a system of M-partial q— difference 
equations. For N — 2, the K— matrix coincides with the one derived in ^U]. One can then construct a XXZ 
spin chain coupled with a quantum mechanical system |22| associated with K^> (it) , in which case the dynamical 
boundary conditions are associated with the AW relations. For general values of N, and any number of sites M, 
it becomes clear that one obtains a XXZ open spin chain coupled with N— type dynamical boundary conditions 
taking if( M+Ar ) (u) in l)18[l. Obviously, the examples above do not exhaust all possibilities. Indeed, starting from 
a different representation for the monodromy matrix satisfying 115|) one gets new quantum integrable models 
which share the same underlying symmetry algebra @. Consequently, from a general point of view this family 
of systems with N degrees of freedom possesses 

- Local integrals of motion derived from the expansion In (i(exp(A)) = ^ n W„A"; 

- Non-local integrals of motion derived from the expansion 1)27(1. 

Applying the first relation in (|18|) , it follows that they are mutually commuting. The existence of such integrable 
structure of the form besides the known local integrals of motion, allows us to say that these models are 
superintegrable. 

• Integrable field theory: Although the construction of Section 3 can not be directly applied to the case 
N — > oo, it is however possible to identify the conserved quantities © in known quantum integrable models. 
For instance, as shown in |23j the sine-Gordon model has a U qo (sl2) symmetry generated by non-local conserved 
charges (up to an overall scalar factor) Q±,Q± satisfying the defining relations ©. In this model, the particle 
spectrum consists of soliton/anti-soliton and breathers. Usually, each particle is associated with an asymptotic 
scattering state \9; m) with rapidity 6 and topological charge 2m where —j<m< j. In this representation, the 
charges act as 1(10(1 setting v = exp((2//3 2 — 1)9) with coupling constant 1 . On a general asymtotic p— particle 
states, they act as A^(Q ± ), A^(Q ± ) 23 . Following [Ej and the results of Section 2, the sine-Gordon model 
or its boundary version (for cq — 1) 24 possess N — > oo non-local conserved charges of form |J3J| in involution. 
The corresponding spectral problem is, in general, rather complicated. However, acting on a p— particle state the 
spectral problem associated with © truncates to the one associated with the conserved charges T^l i , n = 1 , . . . , p. 
In particular, for a single 1— particle state the corresponding eigenstates are expressed in terms of Askey- Wilson 
g-orthogonal polynomials. 

• The limit q — ► 1 and the Onsager algebra: For q = 1 and cq = 1, the structure @ becomes isomorphic to 
the Dolan-Grady one @|, with a trivial coproduct. In the homogeneous gradation of 10} , it is easy to show that the 
loop algebra sl 2 occurs explicitly as A -» E + +£_ and A* -» tE + +t~ 1 E-, t e C, with [E+,EJ\ = 2H, [H, E±] = 
±E±. In terms of the Onsager algebra generators A m , G n , n e Z one has A m — 2t m S+ + 2t~ m S~ which 
gives the identification A — > Aq/2 and A* — * A\j2. It follows that the AW relations (|llfl in this limit become 

^[A ,[A ,A 1 ]]=A 1 -A_ 1 and ±[A 1 ,[A 1 ,A Q ]]=A Q -A 2 . (40) 

Higher order relations are similarly generated. Also, in this limit it is easy to show that the mutually commuting 
quantities @ coincide with the well-known ones As we mentioned in the previous Section, for N = 1 the 
eigenfunctions for l|28|l are related with the Askey- Wilson q— orthogonal polynomials 1)34(1. In the limit q = 1, these 
eigenfunctions reduce to their "classical" analogue, i.e. Wilson, Jacobi, Racah,... polynomials. In this context, it 
is probably not surprising that the Jacobi polynomials arise [57] in the study of the spectrum of the chiral Potts 
model 28 at the superintegrable point [2SI ■ Finally, let us point out that the integrable models discussed above 
(XXZ, sine-Gordon) at this special point q = 1 enjoy the Onsager symmetry algebra. 

6 Concluding remarks 

In this paper, we have proposed a generalization of the Dolan-Grady construction 0] based on the tridiagonal 
algebraic relations l|2|l. This construction was possible due to the realization of the TD pair in terms of tensor 
products of Ugi/z^sfa), which plays a crucial role in the analysis. It was argued that the corresponding finite set 
of mutually commuting quantities can be either derived by imposing the algebraic structure P|. or generated 
from solutions ((16(1 of the reflection equation. Finally, the tridiagonal algebraic symmetry has been exhibited in 
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various examples of quantum integrable models with finite N (or infinite) degrees of freedom, which are known 
to possess already N local conserved charges. Having identified the fundamental generators A^ N \ and the 

corresponding N (or infinite) conserved quantities to non-local charges (JHJ) in involution, we conclude that these 
models are superintegrable. 

The next step is the explicit construction of the corresponding q— Onsager symmetry algebra in the spirit 
of 30] |2H] that we leave for further investigation. Let us mention that beyond the explicit construction of the 
mutually commuting quantities Z^n-i with 1 < n < N, a new (finite) set of algebraic relations generalizing 
the Askey- Wilson ones can be derived directly from i|13|l . Given N, these relations are responsible of the 
truncation of the integrable hierarchy i.e. any quantity I^n—i f° r n > N + 1 can be expressed in terms of l^n—i 
with n < N. For N = 2, it is not difficult to obtain such algebraic relations satisfied by A^ 2 -*, A*'' 2 '' using the 
explicit form of the off-diagonal entries of K^ 2 \u). We intend to discuss these relations, as well as the related 
spectral problem which arises in the classification of multivariable q— orthogonal polynomials separately. 

To conclude, as the Dolan-Grady relations arise explicitly in various problems [23 GH1 1321 we believe our 
construction will provide a new tool in order to derive exact results in massive quantum integrable models. In 
this direction, it would be interesting to find the extension of the Bazhanov-Lukyanov-Zamolodchikov program 

m 
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